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Abstract. The Tracy-Widom distribution that has been much studied in 
recent years can be thought of as an extreme value distribution. We discuss 
interpolation between the classical extreme value distribution exp(— exp(— x)), 
the Gumbel distribution and the Tracy-Widom distribution. There is a fam- 
ily of determinantal processes whose edge behaviour interpolates between a 
Poisson process with density exp(— x) and the Airy kernel point process. This 
process can be obtained as a scaling limit of a grand canonical version of a 
random matrix model introduced by Moshe, Neuberger and Shapiro. We also 
consider the deformed GUE ensemble, M = Mq + y/2SV , with Mq diagobal 
with independent elements and V from GUE. Here we do not see a transition 
from Tracy-Widom to Gumbel, but rather a transition from Tracy-Widom to 
Gaussian. 



1. Introduction and results 

1.1. Introduction. In the random matrix litterature tliere has been alot of discus- 
sion about the transition from Poissonian to random matrix eigenvalue statistics, see 
for example 0)01 COli [22] i [HI' One motivation comes from disordered systems, 
and another from quantum chaos where Poissonian statistics is expected to describe 
the eigenvalue statistics of classically integrable systems (Berry- Tabor conjecture), 
and random matrix statistics should describe eigenvalue statistics of systems whose 
classical dynamics is fully chaotic (Bohigas-Gianonni-Schmidt conjecture). Hence 
it has been natural to look at transitions between Poissonian and random matrix 
statistics. In general there could be many ways to go between different ensembles, 
but it is nevertheless interesting to find natural interpolating ensembles and investi- 
gate their properties. Mathematically it is easiest to consider Hermitian (Unitary) 
ensembles. Previous papers on the problem have been mainly concerned with the 
transition statistics for eigenvalues in the bulk of the spectrum. In the bulk we 
should see a transition from a Poisson process to a sine kernel determinantal point 
process, and for the nearest neigbour spacing statistics we should see a transition 
from the exponential distribution to the Gaudin distribution. 

In this paper we will discuss the edge behaviour of the eigenvalues. In a finte 
random matrix ensemble we look at a scaling limit around the largest eigenvalue 
instead of in the bulk of the spectrum. If we take a diagonal matrix with indepen- 
dent Gaussian entries, the largest eigenvalue will, as the size of the matrix grows, 
fluctuate according to the Gumbel extreme value distribution. If we take a full 
Gaussian matrix from GUE, then the largest eigenvalue will fluctuate according 
to the Tracy-Widom distribution. Can we find interesting distributions that in- 
terpolate between Gumbel and Tracy-Widom? Is there a family of determinantal 
processes that interpolates? Should we typically expect to see a transition from 
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Tracy- Widom to Gumbcl. To shed some light on these question we will discuss two 
interpolating random matrix ensembles on Hermitian matrices. 

1.2. The ensembles. (1) Deformed GUE. Let diag(2/i, . . . tUn) denote the diago- 
nal matrix with elements yi, . . . , j/at. We take . . . , j/at to be independent Gauss- 
ian, N (0, 1/2) say, random variables. Let V be an independent N x N GUE matrix 
with density exp[— Tr y^JdT^. Consider the random matrix 

(1.1) M = dmgiyi,...,yN) + V2SV, 

where S* > is a parameter. When 5 = we have a diagonal matrix with indepen- 
dent entries and when S — > oo, the matrix M/^/S approaches a GUE matrix. 
(2) MNS-model. This model was introduced by Moshe, Neuberger and Shapiro in 
|21| . and we will call it the MNS-model. Let _ff be a Hermitian matrix and U a 
fixed unitary matrix. A probability measure on the space of Hermitian matrices is 
defined by 

PN,u[H)dH = -^e-Trff=g-6TY([C/,//][a,//]*)^^^ 

where 5 > 0, dH is Lebesgue measure on the Euclidean space of Hermitian matrices, 
[U, H] = UH — HU and the star denotes Hermitian conjugate. The weight is 
maximal when [U, H] = so that U and H can be simultaneously diagonalized. 
The unitary matrix selects a preferred basis. We can get a unitarily invariant 
measure by averaging over the unitary group with respect to the Haar measure, so 
that we consider a random preferred basis. We obtain the probability measure 

(1.2) P^[H)dH=^e-^'"" { I e^^T, ([c/,//][[/,//]-)^f;\ 

\Ju(N) J 

The integral over the unitary group can be evaluated using the Harish-Chandra 
or Itzykson/Zuber integral and this makes it possible to compute the eigenvalue 
measure induced by \1.2\ . This gives 

(1.3) pN{x)d''x = 4r det fe-(''+i/2)K^+-?)+2;.....) ^ 

Zpf V / l<i,j<N 

where xi,...,xn are the eigenvalues of H. Actually we will consider a grand 
canonical version of the model, see below, as was also done in There is a 

generalization of the MNS-model to Laguerre/Chiral type ensembles that we will 
not discuss here, see |7]. 

Both of the models above have interpretations in terms of non-interesecting 
paths. 

(1) Deformed GUE. Consider N standard Brownian motions on the real line, we 
think of them as particles, started at j/i, . . . , at time and conditioned never to 
intersect. Let xi, . . . ,xn be the positions of the particles at time S. The probability 
distribution of xi, . . . ,xn is the same as the eigenvalue distribution of M in Hl.l|l . 
sec for example |12j . 

(2) MNS-model. Consider N standard Brownian motions on the real line started at 
xi, . . . , xn at time 0, conditioned to come back to xi, . . . , xn at time t and without 
having had any collisions during this time. Put an initial density Jli3=i ^n the 
points xi , . . . , xjv • By a theorem of Karlin and McGregor, JJj , we get a probability 
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density 

(1.4) ^dctfe-^(=^-+^')-^(^'^^^)'') d^'x, 

Zn V / l<iJ<N 

on the Xj :s. This is the same as ()1.3|l if we take b ~ l/2t. We can think of this 
as a model of non-intersecting paths on a cylinder. As stated above we will be 
interested in the transition at the edge of the spectrmn. The transition in the bulk 
of the spectrum in the the MNS-model occurs when b/N'^ ~ c or i ~ l/2cN'^, 
c > a constant, as oo. (The bulk transition in deformed GUE occurs 

when S ~ C/N'^.) It is remarked, but not discussed further, in |2] that when the 
bulk transition occurs, the behaviour at the edge is still like that of independent 
eigenvalues. Below, we will see that there is a transition at the edge when b/N^^^ ~ 
c as iV — !■ oo. 

1.3. The Gumbel and Tracy- Widom distributions. Consider N independent 
random variables Xi, . . . , Xn with distribution N(0, 1/2). Then it is well known 
that, .2ii|, 

max(Ari, .... Xn) — 



(1.5) I 
as iV — > oo, where 



(1.6) flAT = VlogiV- 



< X Fg{x) = e 
log(47r log A^) 

4VT^iAr 

1 



2VlogA^ 

The distribution function Fq is often called the Gumbel distribution. If we think 
of Xi, . . . ,X]\[ as a point process on the real line with A^ points and we take the 
appropriate scaling limit around the rightmost point we get a Poisson process on 
R with density e^^. Its correlation functions are 

k 

(1.8) pk{xi,...,Xk) = Y[e-''\ 

i=i 

fc > 1. 

The Tracy- Widom distribution, Ftw, is defined by the Fredholm determinant 

(1-9) Ftw{x) = det(/ - KMry)L^(x.,co), 

where 

/>oo 

(1.10) KMry{x,y)^ Ai (x + A)Ai (y + A)dA, 

Jo 

is the Airy kernel, '30'. This distribution occurs in several different places and has 
been much studied in recent years, see jl4j . |31j for reviews. 

If we have a point process on R then its correlation functions, pk{xi, . . . , Xk), are 
characterized by 

(1.11) E[n(i+'/'(^,))]-E M / n'^(^j)'°fc(^i'---'^fc)^'^' 

j fe=o ■ j=l 

for any measurable bounded function on R with compact support. Here the product 
in the right hand side is over all particles in the process. 
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A point process on K. is called determinantal if all its correlation functions, 
Pk{xi, . . . ,Xk), k > 1, exist and are given by 

(1.12) pkixi, . . . ,Xk) = det{K{xi,Xj))i<ij<k 

for some function : ^ IR, the correlation kernel. A Poisson process on K with 
density p{x) can be viewed as a, somewhat degenerate, determinantal process with 
correlation kernel, 

(1.13) KcAx,y) 

1.4. The interpolating process. Before we discuss the asymptotics of the MNS- 
model and the deformed GUE model we will consider a determinantal process which 
interpolates between the Poisson process with density e~^ and the Airy kernel point 
process, i.e. the determinantal process that has kernel (|1.1UI) . We will see later that 
this process can be obtained as a scaling limit of the (grand canonical) MNS- model, 
and we will call it the interpolating process. Also, we will see that we do not have a 
transition between the Tracy- Widom and the Gumbel distribtions in the deformed 
GUE ensemble. Rather we will see a transition from Tracy- Widom to Gaussian. 
This will be discussed further below. 
Define 

/■oo a\ 

(1.14) M^{x,y)^ ——^Ai{x + X}Ai{y + X)d\. 

That the integral is convergent follows for example from H2.ll) and the Cauchy- 
Schwarz inequality. 

Proposition 1.1. The kernel defines a trace class operator in L^(a,oo) for 
any real a, and there is a determinantal process, the interpolating process with 
correlation kernel M„. 

That the kernel interpolates between the correlation kernels for the Poisson 
process with density and the Airy kernel point process is seen in the next 
theorem. The theorem will be proved in section |21 

Theorem 1.2. We have the following scaling limits 

(1.15) lim -M^{-~^\og{Ana''),-^^\og{ATTa'')) = K,,t{u,v) 
a^o+ a a la a Za 

and 

(1.16) Ma{u,v) ^ KAiry{u,v). 

a — ^oo 

It is not hard to see that Ma{x,x)dx < oo for any t, see ()2.2|l below, and 
hence the interpolating process has a last particle almost surely. The distribution 
function for this last particle will interpolate between the Gumbel and the 
Tracy- Widom distributions. 

Theorem 1.3. The distribution function for the last particle in the interpolating 

process is 

(1.17) 

i-iy 



oo 



Fait) = y j / det{Ma{xi,Xj))i<i_j<nd'^X = det{I - Ma)L^(t,oc)- 

^0 -^(i.-)" 
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(1.18) lini F„(i - -L log(47ra3)) ^ Fg{0 
and 

(1.19) lim F„(0 =Ftw(0- 
We postpone the proof to section |2J 

There is a different way of obtaining the distribution Fa that is given in the next 
proposition, which will be proved in section [3 The construction in the theorem 
will not give us the whole interpolating process though. 

Proposition 1.4. Let xi > X2 > ■ ■ ■ be a realization of the Airy kernel point 
process. Let 1/1,1/2, ■■ ■ be independent random variables with common distribution 
function 

(1-20) G„(x) = 

which are also independent of {xi}. Define a new point process by Zj — Xj + yj, 
J > 1. Then, 

(1.21) P[maxzj < = -Fa(C)- 

Note that the point process {zj} in the theorem is not the interpolating process, 
it is only the last particle distribution that is the same. Processes with Poissonian 
edge behaviour constructed in a similar way have recently been studied in j25| . 

1.5. Grand canonical determinantal processes. The probability measure (|1.4|l 
on does not define a finite determinantal point process on M. To get a determi- 
nantal point process we have to consider a grand canonical ensemble with varying 
A'', see for example 0], jJS] for related constructions. Let us first consider a general 
model with the same structure. 

Let X be a complete separable metric space with a reference measure /i. Assume 
that j > 0, is an orthonormal family of complex- valued functions in L'^{X, /i). 
Also, let a„ > be a sequence such that J2'^=o ^ 

00 

(1.22) (j){x,y) = ^an■^l;n{x)■^pn{y)■ 

n=a 

The function (j){x, y) is well-defined in L^(/i x fi) and 4>{x, x) is well-defined in L^{fi). 
We can define a probability measure on X" by 

(1.23) pN{x)d^fi{x) = det{(j){xi,Xj))i<i^j<Nd'^ fi{x), 
where 

(1.24) Zn = det{4){xi,Xj))i<i,j<Nd'^ n{x). 

Here we assume that pn{x) > and > 0. We construct a grand canonical point 
process, compare |2j, p. 123 , by letting 

(1.25) QN 



Nl Z{X) 
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be the probability of seeing exactly N particles, and pN {x)d^ fj.{x) be the probability 
measure for finding particles at xi, . . . , xjv given that there are exactly N particles. 
Here Z{X) is a normalization constant (grand canonical partition function), 

oo , ^ 

(1-26) Z{X) = m^^' 

where Zq — I. If g is a function in L°° with bounded support, then 

oo „ JV 

(1.27) E[n(l+.g(.^,))]= Ef^ / n5(=^^^)det(0(a^.,^,))i<M<ivd^M^), 

j N=0 ^ j=l 

where the product in the left hand side is over all particles in the process. The 
next theorem, that will be proved in section O says that this construction leads to 
a determinantal process. 

Theorem 1.5. The grand canonical point process defined above is a determinantal 
process on X with correlation kernel 

oo > 

(1.28) i^,(a;,y) = ^^^V42:)^n(y)- 

We call this type of process a grand canonical determinantal process. 

1.6. The MNS-model. In the MNS-model we have X = K and ^ is the Lebesgue 
measure, in the above construction. We take 

(1.29) ^{x, y) = 0t(x, y) = ^e''^-'W)/2H--vf /2t 



Then the probability measure (|1.4|l is exactly the measure p.23|l . That we have an 
expansion of the form (|1.22|) follows from the next lemma, which is just a way of 
writing Mehler's formula. We will give the details in the beginning of section^ 



Lemma 1.6. Set f3g = \[^- Th 



en 



— exp(-i(x2 +2/2) - TT^ix - y)'^) 



(1.30) = 5^/3,g"+i/2/j„(;3,)/i„(/3,y)exp(-^(a;2 +2/2))^ 

n=0 

where hn{x), n > 0, are the normalized Hermite polynomials. 
If we make the identification l/2t = q/{l — qf' and define 

(1.31) M^) = y^/i«(/3gx)e"'''"'/2. 



then the ipm n > I are orthonormal and (j) = 4>t can be expanded as in p.22|l with 
a„ = g"+i/2 Theorem II . 51 then gives the next theorem. 

Theorem 1.7. The grand canonical MNS-model coming from ^.!^) or il.4\ ) is a 

determinantal point process on R with correlation kernel 

(1.32) Kx{x,y) = ^ +\ n+i/2 M^)My), 

n— 
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where the three parameters are related by 
(1.33) 6=1 = 



2t (l-g)2- 
Write 

(1.34) q = e-f' 

and fix a number N > 0. Note that n ^ oo corresponds to i ^ oo and /i 0+ to 
t 0+. If we choose 

(1.35) A = e^^-1, 

then J^ K\(x, x)dx ~ N, so the expected number of particles in the process is 
approximately N. 

The next proposition shows that the kernel K\ interpolates between a point 
process defined by N independent Gaussian random variables and GUE as we 
should expect. We postpone the proof to section 0] 

Proposition 1.8. If we choose q as in \1.S4}) and A as in 1^1.35}) . then 

(1.36) hm K^{x, y) = V /i„(x)/i„(y)e-(-'+^')/2 = KcuEiN) 
uniformly for x^y in a compact set, and 

{-^e~^ if X ~/~ V 
f 
tfx^y 

pointwise. 

As mentioned above the bulk transition occurs when fi ^ 1 /cN. This is the limit 
that was studied and discussed in ,21; . 

Theorem 1.9. Let ji = 1/cN, with c > fixed, and let A be given by 1^1. 35\) . In 

this case X is a constant A = e^/^ — 1. The following limit holds, 

(1.38) hm ( ^) = Ux,y) = H 'S^p^phdu 

uniformly for x,y in a compact set. 

The theorem will be proved in sectional 

Thus in this transition region in the bulk of the point pocess we will have a 
determinantal process with correlation kernel L^- Suitable scaling limits will give 
the sine kernel as c — > 0+ and a uniform Poisson process as c — > oo 

In j21| only the following approximate expression 

n Qn\ T f \^ ^'^ sin7r(x - y) 

^'■^^^ '--^"'^)"Y sinh.M^-,)/2 

is given, valid when c is small. At the end of section^ we will sketch an argument 
leading to this approximate formula without discussing the error. 

As briefly mentioned in [21], but not really discussed, when fj, — 1/cN, the 
behaviour at the edge is still like that of independent particles, i.e. we get a Poisson 
process with density e~^ . More precisely we have the following theorem, which will 
be proved in section^ 



8 



K. JOHANSSON 



Theorem 1.10. Let fj, = 1/cN, c > fixed and A — e^/"^ — 1 as in the previous 
theorem. Set 

/— - log(47rlogiV/A^c^) 



an. 



2„2\ 

dbN = (2VTogiV)^i as m (T^. Then, 



(1.40) lim hNKx{aN{c) +bN£.,aN{c) + bwil) = Kext{^,r]): 

N^oo 

pointwise. 

To get an intermediate process at the edge we have to pick a larger /i. In fact 
the intermediate process will be exactly the interpolating process with kernel Ma 
discussed above. The next theorem will be proved in section^ 

Theorem 1.11. Choose fi = a/N^/^ , A — e"-^^'^ ~ 1, where a > is fixed. Then, 



(1.41) lim ^-^i^A 



uniformly for £^,rj in a compact set. 

Hence, in the grand canonical version of the MNS-model, we can see a transition 
between Gumbel statistics and Tracy- Widom statistics for the largest eigenvalue. 

1.7. The deformed GUE model. We turn now to the deformed GUE model 
(|1.1|) . The bulk transition in this and related models has been discussed for example 
in 0, jJUj, [25 and we will not discuss it here. It occurs for S ~ c/N"^, which is 
the same as for the MNS-model. When we look at the edge, the behaviour of the 
deformed GUE will be different than that of the MNS-model. We will not see a 
transition between Tracy- Widom and Gumbel. If we choose S = a^/N^^^ we will 
see a change at the edge behaviour as we vary a, but the transition will be from 
Tracy- Widom as a ^ cxd to Gaussian as a 0+. Informally we can interpret 
this as follows. The eigenvalue distribution is approximately a semicircle and with 
yi, . . . , yN fixed we would see Tracy- Widom fluctuations. However, the fluctuations 
of yi,. .. ,yN causes the semicircle to fluctuate, that is the position of the edge 
fluctuates like a Gaussian. We can think of the semicircle as fluctuating basically 
like jf J2iLi Vi' a Gaussian. The effect is that the largest eigenvalue will 

fluctuate like a Tracy- Widom random variable plus an independent Gaussian. There 
is some similarity between this problem and the random growth model with random 
parameters studied in 

Theorem 1.12. Let dfj,(t) be a probability measure on M satisfying J td^(t) = 
0, jt^djiit) — and / |t|^d/z(<) < oo. Let yi,...,j/jv be independent random 
variables with distribution dfi(t) and consider the random N x N matrix 



(1.42) M ^ diag{yi,...,yN) + V2SV, 

where S — c? jN"^!"^ and V is an independent GUE matrix with density Z^'^ exp(— TrV^ 
Let X^Jx be the largest eigenvalue of M . There is a number R{N) ~ 2aN^^^, given 
by i5.11\} below, which depends on rf/i, a and N , so that 

A^falc - RiN) 



(1.43) lim 

AT— >oo 



a/VN 



< t 



'[X + Y < t], 
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where X and Y are independent, X has the Tracy- Widow, distribution and Y has 
distribution N{0,a^/a^). 

If we want to compare with proposition 11.41 we can let > X2 > ■ ■ ■ be a 
realization of the Airy kernel point process and y be an independent random variable 
with distribution iV(0, jo?). Set Zj = Xj + y. Then maxj>i xj = xi + y, will be 
distributed according to the right hand side of (|1.43() . 

Remark 1.13. Another model for the transition between independent eigenvalues 
and GUE random matrix eigenvalues is a band Hermitian matrix with Gaussian 
elements. Let ma be independent N(0,l/2), 1 < z < and Remy, Immy, 
1 < * < j < A^i and j ~ i < b, he independent N(0,l/4), for some given b, 
1 < b < N . Set mij — for 1 < i < j < N , and j — i > b and ruij = ruji. 
Then M — {mij)i<cij<N is a diagonal matrix when b — 1 and a GUE matrix when 
b = N. It is conjectured, see for example |TS3 and references therein, that the local 
bulk statistics, in the limit N —^ 00, shows a transition from a Poissonian to a 
determinantal sine-kernel point process when b ~ cN^^"^, < c < 00. When is 
there a transition at the edge? Based on the results above one might guess that 
the edge transition takes place for a larger 6. Do we see a transition from Gumbel 
to Tracy- Widom or is there something else happening in between? It is not easy to 
approach these problems. Since we are dealing with the edge and not the bulk it 
could be that the method of moments, used with great success in for Wigner 
matrices, is useful here also. 

Remark 1.14. The comparison of the Tracy- Widom distribution with the Gumbel 
distribution suggests that we are thinking of the Tracy- Widom distribution as a kind 
of extreme value distribution. One way to motivate this is as follows. Let w(i,j), 
€ Z^, be i.i.d. geometric random variables, and let ttI!^\ k = 1, . . . , {^^), be 
all up/right paths from (1, 1) to (A^, A^). Set 



For A^ large each Xf. ' is approximately normal. Clearly, the X^. ' are not in- 
dependent. The random variable G(N ^ N) = maxj, A:*;^' is thus a maximum over 
dependent random variables each of which is approximately normal. We know, , 
that G{N,N), appropriately rescaled converges to the Tracy- Widom distribution, 
which thus arises as an extreme value distribution for certain dependent random 
variables. We are not aware of any last-passage percolation problems that would 
interpolate between Tracy- Widom and Gumbel. 

In measures on partitions both the Gumbel and the Tracy- Widom distribution 
appear, |32| . Are there any natural measures on partitions that interpolate in the 
way that the MNS-model does? 



In this section we will give the proofs of the results for the interpolating de- 
terminantal process with correlation kernel Ma. A basic identity that is useful is 




2. The interpolating model 
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for a > and all x,y, see for example [SH] . 

Proof. (Proposition ^jj. We first prove that nia defined by H1.14|l is a trace class 
operator on L^{a, b) for — oo < a < b < oo. Note that M is symmetric and 



Zi^jMa{Xi, Xj) — / 3I3rXT 



2i Ai + A) 



i=l 



for any complex numbers zi, . . . , zjv and all 0;^, so M{x, y) is a Hermitian positive 
definite function. Hence by [JHIj it suffices to show that 



(2.2) 



f 

J a 



Ma{x, x)dx < oo. 



It then follows that defines a trace class operator on L'^{a,b) with Tr Ala = 
Ma{x,x)dx. The inequality 12.2|l follows from the estimate 



Ma{x, x) 



oo 



-Ai {x + A)2dA 



(2.3) 



< 



e"^Ai {x + \fd\ 



/47ra 



-a3:+a^/12 



by JTTll . 

If we can show that < A/^ < /, it follows that there is a determinantal process 
with correlation kernel Mq,, see l^^. Let / be a continuous function on the real 
line with compact support. Then, 



oo /"OO 



Ma{x, y)f{x)f{y)dxdy 



oo J —CO 



oo 



e"-^ + 1 



Ai (a; + A)/(x)dx d\ 



by Fubini's theorem and hence > 0. Fix e > 0, < e < a, and note that 

oqA 



(2.4) 



oQA 



1 



< e 



eA 



for all real A. Thus 



oo goA 



oQA 



< 



Ai{x + X)f{x)dx] dX 



OO 

OO nOQ / pOO 



^Ai(x + A)Ai iy + X)] f{x)f{y)dxdy 



— oo J — oo \J — oo 
oo poo -1 



oo ^ — oo 



/Aire 



^{x-yy/4t-e{x+y)+t^/12 



f{x)f{y)dxdy. 



Since / is continuous and has compact support this last integral 
Since e > can be taken arbitrarily small we obtain Ma < I. 



2 as e ^ 0- 



□ 



Next we turn to the scaling limits of the kernel Ma- 
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Proof. (TheoremOJ. Set /(a) = ^log(47ra3). Then 

\ U V 

M^iu, v) = -M4- - /(a), - - /(a)) 
a a a 

I |-oo gQ(t+/(a)) 



(2.5) = %/47ra / ^ Ai(i + u/a)Ai(f + 

J-oo V47ra3eat + 1 

Using the identity H2.1|l this can be written as M*(u, v) — A — B, where 

^ g(M-t))V4ct^-(«+i')/2+a3/12 

/oo ^at 
^ Ai + u/Q:)Ai (t + v/a)dt. 
-oo V47ra3e"* + 1 

It is clear that ^ — > as a ^ 0+ ii u ^ v and A — > e~" as a ^ 0+ if u = w. Hence, 
we have to show that S ^ as a ^ 0+. We can assume that u > v without loss 
og generality and write 

B = B, + B2 + B3=( + + ]^== Ai(^ )Ai(^ )dy. 



Now X — > Ai (x) is a bounded function and we have the estimates 
(2-6) 

for a; < 0, and 



(2.7) |Ai(x)|< 



for a: > 0, where C is a numerical constant. It follows from these estimates that 



Jo 



Vv 



Clearly, Bi as a 0+. Similarly, 



and hence ,62 ^ as a ^ 0+. Finally, 



^2y 



dy, 



which goes to as a ^ 0+. This proves (i) in the theorem. 

The fact that Ma{u,v) KAiry{u,v) as a — > cx) follows from the estimates 
(|2.t)|) , H2.7|l and the dominated convergence theorem. □ 

It follows from the estimate (|2.2|l that the interpolating process has a last particle 
almost surely. Its distribution function is given by theorem ll .31 which we now prove. 
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Proof, f Theorem I I.ISII . It follows from Hadamard's inequality and l|2.6|l that 

/ det{Ma{xi,Xj))d"-x <^— ( Ma{x,x)dx] < oo. 

Hence, the first equality in (|1.17() holds, see for example ^Hl- The second inequality 
follows since Ala is trace class on L?{t, oo) for any t and TrAf^ = AIa{t,t)dt, 

m 

Next, wc turn to the proof of ((T~TF|) . Write /(a) = ^ log(47ra^). Then 
(2.8) /(a)) = ^ / det(M:(x„a;,))d"a:, 



where Ma is as in H2.5(l . It follows from the estimate (|2.3(l that 

M:{x,x) < ^^e-"+"^'(")+"'/12 = g-.+a-Vl2^ 

Hence, by Hadamard's inequality, 

and it follows from H1.15|l . (|2.8(l and the dominated convergence theorem that 

n=0 

which proves (|1.18l) . 

If we use the estimate (|2.4|) with e = 1 and (|2.1(l we see that for a > 1 

(2.9) Maix,x) < e-"". 

Hence, for a large, we have 

det{Ma{xi,Xj))i<ij<n < e"^3=i"=j. 

Consequently, we can use (I1.16f) , p.l7|l and the dominated convergence theorem to 
conclude that 11.19|l holds. □ 

We will now prove proposition II . 41 which gives an alternative representation of 
the Fa -distribution. 

Proof. (Proposition ■ We have 

oo oo 

P[maxz, < - E [n(l - X(C.oo) (^,))] = E [11(1 " Xi^M^j + Vj))] 
= 



n yj j^ - - y))dGa{y)^ 



where Ea, denotes expectation with respect to the Airy kernel point process. Here 
we have used the fact that the j/^ is are independent with distribution Ga and that 
they are independent of the Airy kernel point process. The last equality then 
follows from Fubini's theorem. We have also used the fact that the Ga-distribution 
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is symmetric to replace y with — y. Note that X{i.oo)ixj — y) = if and only if 
y > Xj ~ and thus the last expression can be written 



Y[{1 - G^{x, - 0) 



OO 



since we have the Airy kernel point process with correlation kernel (|1.1(J|) . We now 
make the shift Xj ~> + ^ and manipulate the expressions as follows 

V^-p- / TTG„(xi)det( / Ai{xi + t)Ai{x^ +t)dt)d''x 



fe=0 

OO 

E 

k=0 



i=l 



d''x]TGQ(Xi) / d'^t det{Ai {xi + ti) Ai (xj +U)) 



£ I '^'^ / '^'^ ( n Go.{x^)Ai (x, + to ) det(Ai 



i (Xj + tj)): 



where we have used the fact that the determinant is unchanged under transposition. 
Now, this last expression can be written 



det 



GQ(x)Ai {x + ti)Ai {x + tj)dx d'^t 



det(/-M«)i2(^,o„) ^FaiO- 



□ 



3. The grand canonical point process 

In this section we will show that the grand canonical point process defined in 
section fusing (|1.23() and H1.25|l is a dcterminantal process with correlation kernel 
given by (|1.28() . The proof is based on the identity H1.27|l . 

Proof. (Theorem • We want to prove that 

E [11(1 +5(2^,))] l[gixj)demx.,Xj))i<,,,<Nd''fi{x) 

j ^ ' N=0 ■ j=l 

(3.1) = E 1^ / Y[9ixj)det{Kx{x^,Xj))l<i^J<Nd^^i{x). 

N=0 ■ j = l 

The first equality is just p.25|l and p.27|l . The identity 1)3.111 implies the theorem, 
see e.g. ^6j. To prove H3.1|l we will use some facts on von Koch determinants, see 
e.g. Let (oyOij^o infinite matrix and assume that 



(3.2) 



< OO 



< OO. 



i,J=0 
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Then, 

OO ^ 

(3.3) det(/ + A) = J2- ^ det(a™,„Ji:,-=i 



iv. 

n=0 mGN" 



is well-defined. Furthermore, if we have two such matrices A and B, then 
(3.4) det(/ + A) det(/ + B) = det(/ + A + B + AB). 

Inserting l|1.22() into the left hand side of (|3.1|l we get 

= Z(A) 



iV=0 " ■ ' j=l ' ' ' meN« j=l 

N 



1 °" " f C \ 



mGN™ J 

1 



Z(A) 
where 



dei{I + Dg 



Dg{i,i) ^ ^a]'"^ I {1+ g{x))ipi{x)il)j{x)d^i.{x)a/'^ 



3 

X 



Here we have used the determinatal identity 

(3.5) '^^^^j =^ J ^det{(j)i{xj))det{ipi{xj))d^ fj.{x), 

where all the determinants are of size NxN, see e.g. JB] ■ Clearly, Dq — XuiSij , 
and hence Z{X) ~ det(/ + A), where 

A{i,j) = XuiSij. 

Set 

■Bg(«, j) = Y~p^ J 9ix)->JJt{x)ipj{x)dfi{x)ay'^. 

Note that since g is bounded we have \Dg{i,j)\ < Cfi'J'^a^'^ , \Bg(i,j)\ < CaJ'^aj'^ , 
so the conditions (|3.2() are satisfied. Note also that, 

% + Aai(5y + Bg{i,j) + iABg)^j ^ 5.^ + Dg{i,j). 

Hence, by (|3.4ll 

1 , „ s 1 



det(/ + Dg) = ^^^^^ ^ det(/ + A) det(/ + Bg) 

j<N 



Z{X) ' det(/ + A) 

OO ^ 

= dct(/ + Bg) = ^ — det(Bg(m„mj))i<,j< 



iV! 

N=0 meN" 



°° 1 ^ Aa 1 /■ ^ 

H 1^ H n 1 I T' 1^ / det(V'™.(a;j))det(V'„,(xj)) J|5(2;j)d^^(x) 

Xa 



OO -. / OO 

Ar=0 ■''^ j = l \m=0 



Aa, 



-1pm{Xi)->pm{Xj) d fl{x), 



N=0 

which is the right hand side of H3.1|l . Here we have used the identity (|3.5|l again. □ 
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4. The MNS-model 

In this section we will give the proofs for the results on the MNS-model stated 
in section ^ First we must prove lemma 11.61 which makes it possible to use the 
formalism for a grand canonical determinantal process and obtain theorem 1.6, 
which is the starting point for the asymptotic analysis. 

Proof. fLemma ll.6(l . We will use Mehler's formula, 

4-1) V ^TTl g = n =exp-- j(x^+r) + - ^xy), 

where < q < 1 and Hn are the standard Hermite polynomials. If we use instead 
the normalized Hermite polynomials, 

(4.2) Kix) = 

and rewrite the exponent we obtain 
(4.3) 

Y.q-K{x)KXyy-^'^^^'^^l^ = — =L=exp(--i— i^(.x2+z/) + ^(x-y)^). 

- 9' 2(1 + 9) 1-r 

The change of variables x —* PqX, y — * Pqy and multiplication by I3qq^^^ now gives 

(tonii . □ 

The choice of the parameter A in (|1.35|l with q given by H1.34|l is motivated by 
the fact that the expected number of particles becomes 

~ dx 



f Aq"+ ' f 



fix 



+ Ag"+i/2 1 + Ae- 

= -log(l + A) = TV. 
We turn to the proof of proposition II. 81 
Proof. (Proposition • With q as in (|1.34|) and A as in (|1.35|l we have 

Kx{x,y) = J2 1 + (1 - e-A'^!^ie("+V2-A^)M ^"^'^)^"(^)' 

n— ^ ' 

where ipn is given by p.31|l . We split this into two sums, one from ?i = to — 1, 
called Si, and one from TV to infinity, called S2. Since /3q — > 1 as ^ — > cx) we see 
that Si converges to the right hand side of (|1.36|l . We have to prove that S2 — > 
as /i — + cxD. A useful bound is 

(4.4) K{x)e"^'/'\ < 

for all X, see ^Hl- Hence, 



ni/6' 

N 



which goes to as 00. This proves H1.36|l . 
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To prove H1.37|l consider first the case x — y. We have 

1 < 1 + (e^N _ ^)g-(n+l/2)M < 



Since 



we obtain 



K^{x,y) = i^(g^Ar„i)^-(„+l/2)^ V'nW^«(y) 



n=0 



^ n=0 

^— Ve-("+i/2)A'^„(a;)2 < K^{x,x). 
The formula p.30|l gives 

By letting ^ — > 0+ we get the first part of p.37|l . 
Consider now the second 7^ y. Write 

oo 

y) = (e'^^ - 1) ^ e"("+i/2)''^„(x)V'„(y) 



oo 



(e'^^-l)E 

n=0 

Si + S'2- 



1 + (e^A^ _ l)g-(ri+l/2)M 



By JTHI), 



1 - e-^ 

which — !■ as /i — > 0+. Furthermore 



r!=0 



Zy/n I — e >^ 

by the Cauchy-Schwarz' inequality and HI. 30(1 . We see that the last expression 
as ^ ^ 0+. □ 

We turn now to the proof of the two theorems that concern the asymptotic 
behaviour of the kernel K\ in the regime where we have a transition in the bulk. 

Proof. fTheorem ll.9|l . We will use the following asymptotic formula for the Hcrmite 
polynomials, [H], valid for— l + i5<a;<l — i5, (5>0 fixed, 

(4.5) hn{V2^x)e-"^" = I ' l-_(cos[2n^^(x) - \ arcsina;] + 0{-)), 

where 

F{x) = / — y^dy = — (arccosx — x\/ 1 — x'^). 

J -r 2 
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Set An = A-iei/2c^v. Note that An = \-^+0{l/N) and Pq ^ +0{1/N^/'^) 
as N ^ oo. Write fn{x) = hn{x)e^^ Using the asymptotic formula H4.5|l we 
obtain 



Try 



^ n— 1 



1 - 



Trf3qX 



-1/4 



1 - 



/ TT/3gy 



-1/4 



,1/2 



cos 



X cos 



2nF 



'2nF ( _ i 
\2NV2mJ 2 

2Ny/2^) 2 



arcsm 



2NV2m 



arcsm 



2NV2mJ 



0(1) 



1 ,22, 1 1 

7F ^ A-ie"/^^ + 1 

^ n— 1 



2nF 



7r/3„a; \ 1 

: arcsm 

2NV2^J 2 



2NV2^J 



X cos 



2nF 



f 7rf3qy \ 1 



V2iV\/2cri/ 
Note that the sum 

is bounded in N. We have 



arcsm 



2NV2mJ 



0(1). 



1 ^ 1 1 



2NV2cn 2VrjJV 



o( 



^A^5/2 



and 



2VnN J 



TTxW— +0( 



as A'^ ^ cx), X in a compact set. Hence, 

, i ^ (-1)" 



2N^c ^^2N^c2N^c 



' 2N ^ A-ie"/=^ + 1 \N) 

n— 1 



-1/2 



~\_ ^^^^^ ^ X ^ 2 I ^ /* c 

^ E A-ie»/cA- + 1 (^) - ^V^^ " 2 X 

n— 1 ^ 



cos(7r(a; + y)^ — ) 

OS 7r(a; — ?;)\/t 
A-ie*A + l 



uniformly for y in a compact set as iV — s- cx). If we make the change of variables 
t — V? we obtain Lc{x,y) in H1.38|l . □ 



When we are in the transition region in the bulk, the behaviour at the edge is 
still like that of independent random variables. This is the content of theorem 1 1.101 
which we prove next. 
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Proof, f Theorem 1 1.1 0(1 . We split the kernel Kx as follows 

oo 

1 1 



n=Q 
/Mn 



\n=0 n=Afiv+iy 



Si{x, y) + 5*2(2;, y) + Sz{x, y), 



■>Pnix)-)pn{y) 



where Mjv ~ [(1 — (5)ciVlogiV] with 5 > small. Here An has the same meaning 
as in the proof of theorem II. 91 Note that we have the estimate 



(4.6) 

By TTM . 



1 



1 



•N 



6ivS'i(aAr(c) + bN^, aN{c) + bNrj) 
bN 



exp(-^ - (c^iV^ + 0{Nm - Vf + 0(1)) 



(e-^)^) 



as iV — > CX3. The last identity explains the choice of aAr(c) and bN and we will get 
l)1.40|l if we can prove that 5*2 and S3 both tend to zero as N tends to infinity. 



From the estimates (|4.4|) . (|4.6|l and Pq ^ \/2cN we obtain 



\bNSs,{aN{c) + bN^, aN{c) + bNr])\ < 



E 



^-2n/cN _ 



n=Mjv+l 



ll/6 



as TV 00 provided ^ is sufficiently small. 

We can write V'n(aAr(c) + — i^niV^y) , where 



y 



^{V2^+oij^)){V^-i^V^r' iog(4^ iogiv)+e(2v/i^)-^) 



/2n 



A2( 



For a fixed ^, we see that ?/ > 1 + 5 if iV is sufficiently large and I < n < Mn- We 
can then use the estimate, 



(4.7) 



|/i„(V2nx)e-"" I < — 



^1 _~nF{x) 



for a; > 1 + S, 0, which gives 

\hni'\/2rix)e~ 



/4' 



- „l/4 



where C2 is a numerical constant. If iV is large enough, then y > ^JcNn ^ log 
and we get 

\bNS2{aN{c) + 6Arf , aAr(c) + 6Ar77)| 



< 



1 



2y/N ' VToglV n^/ 
which ^ as ^ 00. 



,^1/4'' 



□ 
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The next result shows that the kernel Mq can be obtained as a scaling limit of 
the kernel K\ 

Proof. (Theorem II. 11 |l . Let An have the same meaning as in the proof of theorem 
HH We have 

An ~ e 

with a negligible error. Also, as iV — > oo, 



Write fn{x) = hn{x)e "^^/^ as above. We have 

1 o(-!-)1f L 



^ n— 

-I X /-Mn-1 A/jv W ^ 



^ / \ fc=-oo k=-MN k=MN + l/ 

where Mn — jN^^"^ \ogN with some fixed 7 > that can be chosen. The as- 
ymptotic contribution will come from S2- Here we use the asymptotic formula 

(4.8) h-kiVm + ^) = ^Ai (. + ^)(1 + 0(^)) 

for |fc| < 77V^/'^logiV and u in a compact set. This formula folows from results in 
[3], see pp. Using this we see that 

j}^o. = i^oo WU^ E g-a(fc-i/2)/^V3^/ i (C + ^^)Ai (r; + ^^) 

-Ai(a; + C)Ai(x + 77)da; = Ma(^,7/). 



, e-"^ + 1 

We still have to prove Ei ^ and E3 — > as iV ^ 00. For Ei we use the 
estimate (|4.4|l . which gives 

^ _C_ 1 1 

fc— — 00 ^ ^ 

which goes to zero as N tends to infinity if we choose 7 large enough. 

For l<x<l + 5we have an asymptotic formula for fn{V2rix) in terms of the 
Airy function, see Pj. Estimates of the Airy function then gives 

(4.9) \fn{V2^x)\ < ^e"™(--i)''^ 
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for some constants c, C, when 1 < x < 1 + S. Since F{x) > c{x — 1)"^/^, we can 
combine this with H4.7|l to see that l|4.9(l holds for all x > 1 . For TV sufficiently large 
this leads to an estimate 



< 



It follows that 

AT 

k=MN+l 

as TV ^ oo if we choose 7 sufficiently large. □ 

We give here a sketch of an argument for the approximate expression (|1.39|l for 
Lc{x,y). Integration by parts gives 

sixvTT{x-y)u 2u e"'/" ^ 
2 [°° simr(x — v)u du 



c Jo TT{x-y) cosh^(u2/2c + a/2)' 

where a = log(l/A). When c is small a w — 1/c. Make the change of variables 
u= 1 + ct. This gives 



-i/c T^ix-y) cosh^(t + ci2/2) 

1 /"°° /, , sin 7r(x - + ct) dt 

[^ + ct) — — , , , 2, 



2 J -00 7^(2; - y) cosh^ t ' 

If we use the addition formula for the sine function and neglect terms containg 
we get 

cos n{x — y)ct sin 7r(x — y) 



Lc{x,y) 



—00 

2 



cosh2t + 1 7r(a; — y) 
7r^(x — ?/)c sin7r(a; — y) ttc sin7r(a; — y) 



2 sinh(7r2c(a; — y)/2) Tr{x — y) 2 sinh(7r^c(a; — ?/)/2) 
We see that as c — > 0+ the kernel Lc approaches the sine kernel. 

5. Largest eigenvalue for deformed GUE 

This section contains the proof of theorem 11.121 Consider N non-intersecting 
standard Brownian motions started at yi, . . . , y^ and conditioned never to intersect. 
If we fix j/i, . . . ,?/jv the particle distribution at time S* is a determinantal process 
with correlation functions 

(5.1) Pm,N{xi, . . . , Xm\ y) = dct (ifjv (xi , a;^ ; y) ) i<i j<m , 

where 

N 



(5.2) K^iu, v; y) = -1^ [ dz f ^s-^.-uf /2S ^ -Q ^f^^^ 

(2TnYS Jr w-z jj^ z - y, 

see e.g. Here 7 is a positively oriented simple closed curve containing j/i, . . . , j/^r 
and r a verical line oriented upwards and not interesecting 7, we place it to the 
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right of 7. The formula (|5.1|l also gives the correlation functions for the eigenvalues 
of the hermitian matrix M given by 

(5.3) M = diag(2/i + ...,2/jv) + \/25y 

with V a standard N x N GUE matrix. If is the largest eigenvalue of M then 
A^x has the same distribution as maxi<j<jv a^j. 

Let ¥y denote the probability measure for yi,y2, ■ ■ ■ and let Vx;y denote the ex- 
pectation with respect to the determinantal process with correlation kernel Kjsf{u, v; y) 
given by H5.2|l . Furthermore, we let P^v — Py ^ Px;y be the product measure. We 
are interested in the distribution function 

N 

F^it) = P[AW <t]= E^[n(l - X(t,^){^j})] 

N 

(5.4) =E,[E,;,[[](l-X(t,oo)(a;j))]]. 

When computing the inner expectation we are considering yi, ■ ■ . ,yN as fixed and 
hence we can work with the correlation functions H5.1|) . 
Fix a number e g (1/7, 1/6) and set 

AAr-{yeM^; \y^\<N',l<i<N}. 

Define a cut-off measure d^N{t) by 

(5.5) dfiN{t) ^ --jj--^^j^X[-N^,N^]{t)dfi{t), 
and the function 



(5.6) Gw(z) 



t 



for z £ C\ [-iV^iV<=]. 

Since by assumption fi has finite 7:th moment it follows that 

which — + as A'' ^ 00. Hence, since the expression in the Ey-expectation in (|5.4(l 
is bounded, we can restrict our attention to A^r and use /^at instead of /i, so we 
regard yi, . . . ,yN as independent random variables with distribution /ijv- Denote 
this probability measure by P^^-* . 

Lemma 5.1. There is a real number Wc — Wcin), which is approximately \J NS = 
aN^/^, such that 

(5-7) G'^iw.) = 

for all sufficiently large N . 
Proof. We have 



[z-tf 



The moment conditions on ^ can be used to see that 
(5.8) \G'Az)-\\<% 
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for real z > 2N'^ say. We see that, for N sufficiently large, G'^ is a decreas- 
ing function in [2N'^,(X)). Furthermore, G'j^{z) ^ as z ^ oo and G'^{2N'^) « 
f/4A^^^ > f/aiV^/'^ Hence, when N is large enough, there is a z = Wc such that 
(|5.7|) holds. From (|5.8|) we see that '-^ l/a^N^^^, which gives the asymptotic 

behaviour. □ 

Set 

« f 

(5.9) My) = ~ ^ 717—712 + NG'^w.), 

j=i ^ 

and 

^ 1 

(5.f0) Vc = wc + SY^ 

f ^ in — 

Furthermore, set 



a 



(5.11) R{N) ^wc + + / "-^dfiNiy) 



and define 



N 



Note that 

(5.13) v, = R{N) + j^sN{y). 
Lemma 5.2. 

(i) There is a constant C such that 

(5.14) EWK(2/)]<C. 

(ii) We have the limit 

(5.15) VarW[sjv(2/)]^aVa2 
as N oo. 

(Hi) The random variable SN{y) converges in distribution to N{0, j o?) as N ^ oo 

Proof. Since the 7:th moment is finite and J td^{t) = we get |/ tdfiN{t)\ < 

(A 

V 

4'"^[rN{yr] - Varf )[r^.(y)] = Varf ^ 



CN-^' < CN-^/'^. The definition of Wc and Gn gives E!)^^[rjv(y)] = and thus 



N[-\d^\w,)~G'^{w,f]<^, 



where the last inequality follows from our moment condition. Since Wc ~ aN^^^ 

we see that the right hand side is bounded. To prove (ii) wc compute 

(5.f6) 
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Writing ?/(l-?//wc)""^ = ?/+y^/wc + - ■ • , wc get / y{\-y /■Wc)^^d^.N{y) = o(l). Sim- 
ilarly, writing 2/2/(1 _y/ti;^)2 = y'^ + 2y'^ /wc + ■ ■ ■ , we get / y^ /{l-y/wcfdfiN{y) = 
+ o{l). Using Wc ~ aN^^^ , the identity H5.16|l now yields (ii). The claim (iii) fol- 
lows immediately from the central limit theorem, since sn is a sum of independent 
random variables. □ 

Define to be the set of all y e An such that |rAr(j/)| < CN'^ and |sjv(j/)| < 
CN^. It follows from lemma lST^ fi'l and (ii) that we can restrict ourselves to y £ B^. 
We will study 

N 

(5.17) F*N{t) = EW[xB„(2/)E.;,[n(l " X(t.oo) (a;,))]] 

instead of F/v (t) given by (|5.4() . Hence it is enough to consider a fixed y in B^ and 
work with 

N 

F*^{t;y)^E,,y[l[{l-Xit,oo){^j))] 
^ f_l)fe r 

(5.18) =y.^Tl^ det{KN{x^,Xj;y))d''x. 

r~, f^- J(t.oo)'' 



k=0 



Our problem is then to investigate the asymptotics of Kiq(u,v]y) for y G -Bat. 
Set 

2 ^ 

/(w) = ^ - ^ + H " 

where we choose the principal branch of the logarithm. The number Vc was defined 
so that 

(5.19) f'{w,) = 0. 
We consider 

Q; OL 

(5.20) u = Vc + i—i= ; V = Vc-^ ri—j=, 



where ^ and rj lie in a compact set. 

To perform a saddle-point argument in the integral (|5.2|l we must specify appro- 
priate contours. Let Ci : [0, oo) B t ^ Wc + t + it, C2 ■ [0, 00) 3 t Wc + 1 — it, 
C3 : [0, 00) B t ^ Wc " t + it and C4 : [0, 00) 3 t ^ Wc — t — it. We want to 
show that we can deform F to Ci — C2 and 7 to C3 — C4 in the contour inte- 
gral 1)5.2(1 . Let Cf be the parts of the contours where we restrict t to [0,A + Wc] 
and let ja ■ [—Wc ~ A,Wc + A] 3 t —A — it, where A > 0. Then 7 can 
be deformed to -f 7^ — in l|5.2(l if A is sufficiently large. From (15.11(1 
we see that R{N) ^ 2aN^/^ and since ^ belongs to a compact set we see from 
(ESI, \sNiy)\ < N" and (|^?^ that u > aN^/'^ for all sufficiently large N. Since 
Re(— (z — uY) < —2uA we see that we can let ^ — > 00 and conclude that the 
contribution from goes to zero. 

Choose F to be a vertical line through Wc- We want to show that the part of this 
line that lies in the upper half plane can be deformed to Ci, and the in the lower 
half plane to C2. Set w ^ Wc + t + iA, < t < A. Then, 

g{t) ^ Re ((w - vf) ^ t^ + 2{wc - v)t -f {wc - vf - A^ . 
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For < t < Ajl, we see that g(t) < -A'^/2 for A large, and when A/2 < t < A we 
have 

g{t) < {wc - v)A + {wc - < -|iVi/6A 

for N and A large. Hence we can deform the upper part of F to Ci in 15.2(1 . The 
deformation to C2 is analogous by symmetry. 

Next, we want to localize the integration to a small neighbourhood of Wc- 

Lemma 5.3. Define gi{t) = Re{f{Ci{t)) — fiwc)) for i = 1,2 and 
gi{t) — —Re{f{Ci{t)) — f{wc)) for i ^ 3,4, t > 0. There is a positive constant 
c so that 



(5.21) g,it) < 



-cN^/H^ ,0 <t <aN^/^/2 
-cN^/H t > aN^/^/2 



for all sufficiently large N . 
Proof. Consider gi{t). We have 

1 1 1 ^ 

9i{t) = 25 [K + - 1^] - ^yc{wc + + 2 E + t-yjf + 1") - Re /(«;,). 

Differentiation gives 

N 



Wc ~'^'c . sr^ Wc~ Uj + 2t 



S f^^{wc~yj+tr+t^ 
and hence ^^(0) = by (|5.1U|) . Hence 



2t^ 



f~{ (^^^ ~ - vj + ty + ^2) ■ 

We know that Wc ~ aTV^/*^ and \yj\ < N". If |i| < aN^/'^/2, we see that there 
is a positive constant c such that g'i(i) < —cN^^-^t^ for < i < aN^/^/2. Since 
51 (0) = we obtain the first part of ^^l^i for i = 1. If t > Q:iVi/s/2, then 
9iit) > —cN^^^ for some postivc constant c, and we obtain the second part of 
(irSTli . 

Consider next 33 (t). Again 53(0) = 33(0) — and we get 

^'^^^ ^ " It {wc-yM^c-y,+tr + t^) 

and we can proceed as above. The functions 32 and 53 are treated analogously. □ 

We also need a local approximation of f{w) in a neighbourhood of Wc- By (|5.19|) 
we have f'(wc) — and we also have 

1 ^ 1 

(5.22) /"K) = ^-E(^z:3^--(^)^ 

by H5.7|l and H5.9|l . Furthermore, 

N 2 

(5.23) /(^n^c) = E 



FROM GUMBEL TO TRACY-WIDOM 



25 



Lemma 5.4. For C e C and |C| < N^^^^, 

(5.24) /K + CaN-'^') = /K) + + 

as N ^ oo, where o(l) is uniform for \(\ < N^^^^ . 
Proof Define i?(C) by 

(5.25) f{w, + C) = /(^«c) + /'(We)C + ^/"(«^c)C' + ^/^'^ K)C' + i?(C)- 

Since Wc ~ aN^/^ and — ?/j | is much, greater tlian 1 for N large enough a Taylor 
expansion gives 

(5.26) |E(C)| < CiV^/^ICI^ 

for Id < 1. It follows from (EH), f??^ and (jF?^ that 
(5.27) 

/(«;, + iCa7V-i/6) = /(^J _ " ,.^(y)c2^-i/3 + ^ ^ ^ + i?(Ca7V-V6). 

Since \rN{y)\ < CN" and |C| < N^'^^ we see that \rN{y)C,'^N-'^/^\ < CN-^/^^. 
Furthermore, by ^1^, \R{CaN~^/^)\ < CN~^/^. We can write 



y ' ^^ + y-g'(---^^-)^^^(i + o(i)) 



as iV — > oo. We see now that (|5.24|l follows from (|5.27() . □ 

It follows from H5.2|l . the definition of / and the change of contours discussed 
above that 



(5.28) = Ctff / / ..^^^e---W..^.-./. 



(27ri)2Q2 Jc\^C2 

Consider z on C3 and ?« on Ci. The other cases are similar. Set z — Wc + (— t + 
it)aN^^/^ , w — Wc + (— r + iT)aN^'^/^ , t,T > 0. It follows from lemma [^31 that 
we can localize the evaluation of H5.28|l to i,T < A^^/^^. By lemma 

/H - fiz) - ^(r + irf -l{-t + ^tf + o{l) 

uniformly for < <, r < iV^/^*. Hence, the contribution to H5.28|l from z on C3 and 
w on Ci is 

dt / dr- 



{2mfa Jo Jo (r + ir) - {-t + it) 

Define 

(5.29) K;,{u, V- y) = e-"~'"+^'^'-^^^/^ Kn{u, v; y). 
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We can just as well use K'^ as Kn- If we argue as above for all parts of the contours 
we get 

Q. (y. 
lim —=Kli{vc+^—;=.,Vc + -q—=) 

(5.30) = [ dz [ du;e"''/3--V3+«z-^-^_^ 

where 7' is the contour given by t + it for i < and — i + it for t > 0, and F' is the 
reflection of 7' in the imaginary axis. 

We have to show that the right hand side of (|5.30fl is really the Airy kernel. 
Let C be the contour given by t + i\t\, t E M.. If we change variables by z = i(, 
w — —iuj, then 7' maps to C and V to ~C and we see that the right hand side of 
H5.30|l becomes 



Jc ^iC^ 



which is the Airy kernel KAiry{£., ??), |15| . 
We have proved 

Lemma 5.5. Define 

(5.31) Kt,{i, ??; y) = -^^n{^c + + 



with K*j^ as in 1^5. 29\) and Vc given by 1^5. 13\) . Th- 



en 

(5.32) KT^{£,,-q]y) ^ KMry{S.,fl) 
uniformly for ^, rj in a compact set and y € . 

To control the convergence of H5.18|l we need some more estimates. 

Lemma 5.6. Fix a constant A. There is a constant C , depending on A, such that 
for ^,ri > ^A and all sufficiently large N we have the estimate 

(5.33) /f^(e,r;;2/)|<Ce-i(l«l^'^+l"l^^^) 
for y e Bn. 

Proof. Deform the contour C3 — C4 to 7 = 71 + 72 + 73 and Ci — C2 to F = 
Fi + r2 + Fa, where 71 : {-00, -aSN^'^/^) 3 t ^ Wc + t + it, j2 ■ {S, 6) 3 t ^ 
Wc - aSN-^/^ + aitN-^/'^, 73 : (Q^(5iV~l/^ 00) 3 t - t + it and Fi, F2, F3 are 
obtained by reflection in the line Rez = Wc- 
Note that for z on 7 and w on F, 

1 7V1/6 
< 



\w — z\ 2aS 
From (Ei2Hl), it^^ and (|OT|) we obtain 

7V1/3 



(5.34) x( f (,'^'=ifM-f(^-))+vN'^''a-^Rc{w-w,)^^^^ 
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On r2 we can use lemma ICTI to get 

Re [/K + "(-J + zt)Af-'/') - /K)] = - 5t^ + o(l). 

This gives 

On we can use lemma to get 



- iVi/e'^ ■ 

The contribution from Fi is analogous. Choosing 6 = for rj > 1 and 5 = 1 
otherwise, we get 

The estimate for the other integral in H5.34() is analogous and the estimate (|5.33|) 
follows. □ 

Define the distribution function by 
(5.35) HN{t;y)^y2^—^ det(XAr(6,C,;2/))d'=e 



k=0 



It follows from lemma 15.51 and lemma 15.61 that 
(5.36) lim HN{t;y)^FTw{t) 

N—*QQ 

uniformly for t in a compact subset and y £ B^. 

If we change from Kn to K'^ in H5.18(l and make the change of variables 
Vc + aS,i/VN, we see that 

F^(t) = Ef) [xBAy)GN{VN{t - vc)/a;y) 

Thus 

F-;,{R{N) + WVa^) = EW[xs„(y)GAr(e - SN{y);y)]. 
We can now use lemma IF?^ and (|5.36() to see that 



lim F;f{R{N) +^a/VN) ^ / Frwi^ - u)hiu)du, 
where ft, (u) = (27rcr^/Q!^)^^/^ exp(— u^q;^/2o'^). This completes the proof of theorem 

EH 

Acknowledgement: I thank Y. Chen for drawing my attention to some papers 
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